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SUMMARY 

Equations are -written, which govern the- motion of a gas -lubricated 
hydrostatic thrust bearing. The nonlinear equations are solved on a high- 
speed digital computer for air as the gas. Systematic investigations are 
made of the parameters involved in order to estimate their effect on 
stability. 

A nonstable case is taken as a standard for comparison, and stability 
is achieved by properly varying any of several geometric and environmental 
parameters. The bearing pad volume and rigidity appear as prime control- 
ling factors of stability. Smaller pad volumes and softer bearings result 
in more stable operation. 

Higher temperatures reduce the weight flow and reservoir pressure 
required for a fixed load and clearance and favor stable operation by 
making a softer bearing. 


INTRODUCTION 

The possibility of using a gas as a bearing lubricant has attracted 
attention for many decades. Application has been confined, however, to 
certain installations where low viscosity was considered to be so Inpor- 
tant that the associated difficulties were not determinative. Such in- 
stallations occur chiefly in instruments such as strain-gage balances, 
gyroscopes, and torque -measuring devices. 

In recent years the second important property of gaseous lubricants, 
excellent thermal stability, has caused a renewed interest in the subject. 
An increasing number of applications is being encountered where extreme 
tenperatures make conventional bearing lubricants impracticable. For 
this reason it becomes desirable to understand better the gas bearing 
and one of itB chief deterrents , instability. 
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The present investigation considers a simple nonrotating thrust 
"bearing consisting of two round flat disks Of finite diameter separated 
by a thin film of air. The air under pressure enters the bearing clear- 
ance space at the center, flows to the periphery, and exhausts to the 
atmosphere. The resulting pressure distribution supports the bearing 
load. 

The nonlinear differential equations governing the motion of the 
bearing when disturbed were solved on a high-speed digital computing 
machine. The solutions provide some information about the effect of 
certain parameters on the stability of the bearing. In addition, these 
solutions provide a means for checking the validity of more simplified 
solutions (e.g., ref. l). 


ANALYSIS 

The simplifications made in order to solve a physical problem are 
generally dictated by the specific information that the solution is re- 
quired to provide. The present solutions are intended to provide some 
information concerning the effect of several parameters on bearing sta- 
bility. Accordingly, an effort is made to select those simplifications 
that might least affect the stability results. * 


Bearing Configuration and Operation 

The geometry and nomenclature of the type of bearing investigated 
are shown in figure 1. Pressurized air comes from the constant-pressure 
reservoir through an orifice into the bearing pad. From the pad the air 
flows through the clearance space to the periphery of the bearing and 
exhausts to the atmosphere. 

If the bearing load is increased slightly, the upper disk in figure 
1 will move downward and decrease the clearance. This tends to decrease 
the weight flow out of the bearing periphery and also, by continuity, 
the weight flow into the bearing pad. Then the pressure drop decreases 
across the orifice so that the pad pressure increases, which tends to 
balance the Increased load. In general, the bearing overshoots the new 
equilibrium clearance and begins to oscillate. This oscillation may lead 
to bearing collapse or may be damped out. The equations governing the 
oscillation are presented in the following section. 


Basic Equations 

The flows from the reservoir to the pad and from the pad to the 
periphery are considered separately and then are related through a con- 
tinuity equation. All symbols are defined in appendix A. 
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The weight flow entering the hearing pad is given by (appendix B) 

/ r-lV/2 1 


/ rziV/ 2 

3600Agp R J./ r J I „ < , 

^ ■ Jba ml '<* V "W ’ * 


r-lW 2 


W 1 = 




■£2£2- ** U - * V ,^2rSi 

r-i / l-y\ ! 

• 5600A8Pr * r l-f r ) , l<f 


In normal bearing operation f will probably be between and 1. 

As the bearing oscillates, however, the orifice may choke (lew pad pres- 
sure) or the flaw may reverse in direction (high pad pressure). Both of 
these situations were encountered in the present investigation. 

The flow outward from the pad radiUB to the bearing radius is computed 
from a modified form of the Navier-Stokes equation. The flow is assumed 
to be purely radial, the pressure is assumed to^be a function of radius 
only, and the inertial forces are considered negligible in comparison with 
the viscous forces. With these assumptions, the following equations are 
obtained as developed in appendix B: 


3600rtgh 

* To 

In — 


(Ps 2 - P l) 


3600tfgh 3 p| / 2 Pg 

— If - 

dRT 0 In -S. \ PR 


(B12) 


*(rl + r a r s + rf) 


- P a ) 


it(r| + r a r s + rj)p T 
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MS 


(B16) 


To preserve continuity, the difference between the weight flows into 
and out of the bearing is equated to the rate of increase of the weight 
in the bearing pad. Considering the change of state in the pad to occur 
adiabatically gives the result (appendix B) 

_ crvv 2 ( f NVr /h + d df . dh\ 

W 1 “ w 2 - 3600jtr s p s ,i •& + &£) (2) 
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where the subscript i denotes conditions at initial equilibrium (time 
t = 0 ). 


Now equation (2) with the substitution of equations (l) and (B12) 
gives one differential equation in f and h. A second differential 
equation is obtained by relating the forces acting on the bearing to the 
weight and acceleration of the bearing: 


W d^h 

12g dt 2 


_ _ dh 

F » It 


( 3 ) 


or, by using equation (B16), 

w a 2 !! *( r B + r a r s + L s.\ - an 

T5i 52 3 - i£) ' F ' “ « 


where F is the bearing loading function and a(dh/dt) is a drag term. 
The value of a is taken to be zero in the present analysis. 

Equations ( 2 ) and (4) provide a system of two nonlinear differential 
equations in h and f. These equations were solved for various condi- 
tions on a high-speed digital computer. 


General Considerations 

In the present analysis the bearing is assumed to be in a state of 
initial equilibrium (denoted by subscript i) under a constant load F d 

(fig. 2(a)). At time t » 0, a change in the loading function occurs 
(usually a step change to a value Fj) which puts the bearing tenporarily 

out of equilibrium. The bearing then begins to hunt for the new equilib- 
rium clearance hf corresponding to the new load ‘Ff. If the bearing 
settles out (fig. 2 (b)) or oscillates with constant amplitude, it is said 
to be stable. If it oscillates with increasing amplitude, it is unstable. 
It is convenient to plot the behavior of the bearing as a dimensionless 
clearance parameter (h - h^)f(h^ - hp) as shown in figure 2 (c). 

As a standard example, an unstable bearing configuration is selected. 
The pertinent bearing parameters are then varied one at a time in order 
to examine the ability of each to stabilize the bearing. The standard 
example has the following characteristics : 
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r a , in 6 

*e> ia 

A, sq in 0.003 

4, in 0.008 



W, lb 32.172 

p a , lb/sg In 14.7 

T, °B 530 

li, n> 500 

h ± , in 0.002 

F f , 3i 550 


With the characteristics selected, a solution proceeds as follows. First, 
certain preliminary computations are made: 

(1) The initial pad pressure p s ^ is found from equation (B16) by 
using the prescribed values of F-j_, p a , r g , and r a . 

(2) The initial weight flow ^2,1 is f oun(3 - from equation (B12) by 
using the prescribed value of hi and the computed value of Ps,i* 

(3) Since v i f ± = w 2,i ft* 1 ® sy 8 ^ 111 is in equilibrium), the required 
reservoir pressure p^ can be found from equation (l). 

(4) The pressure ratio f^ is then obtained from p s ^ and p^. 

( 5 ) Although final equilibrium conditions may never be achieved, the 
values of p S} f, hf, and Wf, which would, be obtained under a constant 
disturbance load Ff, are computed as above, except that p^ instead of 
h is now fixed. Second, with initial values of hi and fi now known, 
differential equations ( 2 ) and (4) are solved numerically to yield the 
bearing motion shown in figure 3. 

In subsequent solutions the same process is carried out with Fi 
and hf the same as in the standard example and with one of the bearing 

parameters changed. In this manner, systematic changes in important 
parameters can be made in an effort to improve the stability character- 
istics over those of the standard example. It is important to notice that 
no viscous damping force is included in the equations although such a 
force will be present in any real system. For this reason, the present 
analysis represents an extreme situation. 
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RESULTS 

The results are presented mostly as plots of the bearing clearance 
parameter (h - hf )/ (hi - hf ) against time as important bearing par ame ters 
are varied. Each plot contains the curve for the standard solution as a 
basis for comparison. The letters in parentheses on the plots are for 
use with table I which gives the results of preliminary calculations and 
the values of inportant parameters for each case. 


Effect of Pinal Bearing Load Ff 

The first parameter investigated was the disturbance load F f . Final 

loads of 501, 550 (standard case), and 600 pounds were applied to a 
bearing having an initial load of 500 pounds in each case. The results 
are shown in figure 4, 

Figure 4(a) shows the variation of the actual clearance h, while 
figure 4(b) shows the variation of the clearance p aram eter 
(h - hf )/(hi - hf). Since hf varied from one example to another, a 
clearance parameter such as this one is useful for comparison purposes. 

The curves of figure 4(b) are in close agreement over the first half 
cycle and deviate increasingly from there on. Since the agreement was 
sufficiently close, only one value of F f (550 lb) is used throughout the 
remainder of the analysis. 


Effect of Bearing Parameters 

Pad depth . - It is generally believed (ref. 2) that reducing pad 
volume reduces the tendency to instability. A convenient way to do this 
is to make the pad depth smaller. In the standard case the pad depth is 
0.008 inch, and the bearing is seriously unstable. Values of d of 
0.004 and 0.002 inch are shown in figure 5 along with the standard case. 
The effect on clearance parameter variation is marked. Stability is 
actually achieved at da 0.002 inch. 

For sufficiently large values of d, the influence of pad depth is 
felt only during a displacement of the bearing, and its value does not 
affect the values of any other parameters at an equilibrium position 
(see table i). In this case, pad depth can be changed without altering 
the design operation of the bearing in any way. 

Orifice area . - Increasing the area of the orifice joining the 
reservoir to the pad decreases the orifice .resistance and reduces the 
reservoir pressure (table i) required to pass the initial weight flow. 


p 


4579 



HACA TN 4095 


7 


(Remember that the initial load, pad radius, and clearance and, hence, 
the initial pad pressure and weight flow are fixed.) Thus, f i => p g i /p R 

is increased as orifice area increases. Also, it should he noted that 
AB’/^h, a ratio comparable to a spring constant, decreases as the orifice 
area increases, which indicates a less rigid bearing (table I). 

Figure 6 shows that the bearing is stabilized in this case by mi 
increase in orifice area of less than 0.002 square inch. The bearing is 
quite stable at an area of 0.005 square inch. 

Throughout the analysis, high values of fi and the corresponding 
softness of the bearing proved conducive to stability. On the other hand, 
high values of f ^ diminish the mavl imm equilibrium load that a bearing 

can carry with a fixed reservoir pressure. This may result in a design 
compromise between range and a degree of safeness with regard to stability. 

Fad radius . - A decrease in the pad radius has a twofold effect on 
the bearing stability. The pad volume is decreased, and the initial 
pressure ratio f ± is increased. Figure 7 indicates the variation of 

fi with pad radius for several values of orifice area. A number of fac- 
tors contribute to the effect on fi. As pad radius decreases, the value 
of pad pressure necessary to support the fixed initial load increases 
slightly as shown in figure 8(a). This variation is determined completely 
by the variation of r s (see eq. (B16)) and is independent of orifice 
area. From equation (B12) the weight flow varies with r s explicitly 
and Implicitly through p B . Figure 8(b) shows this variation of weight 
flow. The rapid increase at large pad radii results in a rapid increase 
in reservoir pressure and, hence, the rapid decrease in fi indicated 
in figure 7. At the smaller orifice areas the effect of the weight flow 
change on the reservoir pressure is more marked. 

Figure 9 shows the effect of decreasing pad radius on stability. As 
would be expected, the effect is favorable, and stable conditions are 
achieved with a pad radius of 3 inches. At a pad radius of 2 inches the 
oscillation is damped very quickly. 

Temperature . - It has been shown that a higher temperature increases 
the load capacity of a bearing at a given weight flow (ref. 3). Or, as 
in the present case, the weight flow required to support a given load is 
reduced (eq. (B12)). The reservoir pressure is consequently reduced, and 
the value of f-j_ is increased. As shown in table I, AF/Ah is decreased, 

which indicates a softer bearing. Figure 10 shows a favorable effect on 
stability. The general improvement resulting from increased temperature 
appears more significant in the case of higher loaded bearings, as is 
indicated in the section entitled "Higher loading." 
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Mass of system . - Figure 11 shows the effect of increasing the com- 
bined mass of hearing and load by a factor of 10. The starting frequency 
is decreased from about 285 to about 98 cycles per second. Thus, the 
frequency varies nearly inversely with the square root of the mass as in 
the case of a linear spring. The first half cycle dip is not as low for 
the larger mass, but the second half cycle peak is higher. The tendency 
toward instability is not appreciably affected in this example. 


Other Effects 

Linear applied load . - The step load used so far in the analysis 
represents an extreme condition never realized in practice. In order to 
gain an appreciation of the effect of this simplification, a number of 
linear load functions were used, each of which distributed the applied 
load over a longer time interval. The results, represented by the two 
cases shown in figure 12, indicate that there is no elimination of in- 
stability in this example. The bearing can be brought down quite close to 
its final equilibrium position if the load is applied very slowly. As 
soon as the load is fully on, however, and the forcing function is constant, 
the bearing begins to oscillate with increasing amplitude. It seems that 
the magnitude of the applied load and the way in which it is applied have 
little effect. 

Isothermal compression . - Throughout the present analysis the air 
compression in the pad has been assumed to occur adiabatically. As this 
may not be true in the real case, a solution of the standard example was 
made in which the pad compression was assumed to occur isothermally, and 
the results are shown in figure 13. The isothermal case (appendix B) is 
the more unstable of the two. 

Higher loading . - In many applications where an air bearing might 
prove useful, the bearing is subject to quite high loading. If the stand- 
ard case is considered with the initial load and step load Increased, the 
instability is heightened as shown in figure 14. The higher load raised 
the pad pressure. Since the initial clearance was the same, the weight 
flow was increased. Therefore, the orifice operated nearer to a choked 
condition (see table i), and h was higher. 

This may mean that the bearing must be designed for the maximum load 
that it will be expected to carry. Figure 15 shows the result of a very 
highly loaded (15,000 lb) bearing having the same outside diameter as the 
standard e xam ple with pad radius, depth, and orifice area chosen to elimi- 
nate instability. The small pad depth indicated in table I is achieved 
practically by making an annular pad instead of a circular one. The ef- 
fective circular pad depth can be made a small fraction of the actual 
annular pad depth. 
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Very often high temperatures are present in the situations where the 
bearing loads are high. This is a help in reducing the weight flows re- 
quired and in obtaining conditions conducive to stability. Figure 16 
illustrates the favorable effect of increased temperature on a highly 
loaded bearing. The respective weight flows and reservoir pressures are 
indicated in table I. The higher temperature reduces the required weight 
flow to less than one-sixth and the reservoir pressure to less than one- 
third the original. The effect on stability is primarily due to the 
increase in fi and consequent softening of the bearing. 


CONCLUDING REMAKES 

In view of the simplifications made for the preceding analysis, ex- 
perimental data are needed for substantiation of the results. Although 
such data are scarce at the present time, this situation may change in 
view of the potential that gas lubrication has for bearing applications 
of the future. 

Certain conclusions can be drawn, however, from the analysis. Gas- 
lubricated bearings can evidently be designed to be stable, even under 
very high loads. Proper choice of pad volume and bearing rigidity is 
necessary. With the help of carefully obtained experimental data, either 
the derived equations or simpler linearized equations can probably be 
used to determine stable ranges of the important parameters. TTl gb tem- 
peratures are useful to the designer in lowering weight flow and reservoir 
pressure and raising the ratio of pad pressure to reservoir pressure. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 31, 1957 
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APPENDIX A 
SYMBOLS 

A effective orifice area, sq in. 

d depth of pad, in. 

F force exerted on hearing by applied load, lb 

F-g force exerted on bearing by fluid, lb 

AF F^ — F^ , lb 

f ratio of pad pressure to reservoir pressure, P s /pp 

f c ji ratio of pad pressure to reservoir pressure -when orifice is choked 

g acceleration due to gravity, ft/sec^ 

h bearing clearance, in. 

Ah h^ - hp, in. 

p static pressure, lb/sq in. abs 

R gas constant, ft2/(sec^) (°R) 

r radius , in . 

T temperature, °R 

t time, sec 

u radial velocity, ft/sec 

V s effective pad. volume, Jtr|(h + d) , cu in. 

W weight of combined bearing and load, lb 

w weight flow, lb/hr 

w^ weight flow into bearing pad, lb/hr 

w 2 weight flow out of bearing pad, lb/hr 
z coordinate in direction of bearing axis. In. 
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a ratio of drag to velocity, lb-sec/ft 
8 time interval, sec 

X ratio of specific heats 

M- viscosity, lb-sec/sq. ft 

p density, Xb/cu in. 

Subscripts : 

a condition at outer edge of bearing 

f condition at final equilibrium, if reached 

i condition at t = 0 

R condition in reservoir 

s condition at outer edge of pad 
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APPENDIX B 


DEVELOPMENT OF EQUATIONS 

The weight flow through an orifice from a reservoir of higher pres- 
sure to one of lower pressure p^ is given "by the equation (e.g., 

ref. 4) 


w = 360QA ■ 





Using the perfect gas law 


P = 


12pRT 

g 


(Bl) 


(B2) 


in equation (Bl) and simplifying give 

3600Ag /P L \ l/r 

' Ph V p h) 


4 £U 

y 2r 


L H 


r-r 

r 


i - 


|l/2 


(B3) 


where Tg is the temperature in the higher pressure reservoir. If 
p s = Pg in the case of the gas -lubricated "bearing. 




3600Ag 
RT = 


Ps 



i fr 


r-r 

r 


1/2 


(B4) 


where the minus sign indicates that the flow is from the pad into the 
reservoir. With the notation f = p s /p equation (B4) can "be written 


i-rW 2 


w, ■ . -^ -- Y^U-f r) 


1 Jr-- A 

y 2 r 


RT 


(B5) 


8 


* 


4579 



4579 


NACA TN 4095 


13 


s> 


4 . 


If 


p s = p R , tut the orifice is not choked, 


equation (B3) "becomes 


w. 


5600Ag 




(B6) 


If the orifice is choked, the constant weight flow is given by 


w 


1 


5600Ag 


i 

r 

P R f ch 



ni \ 1/2 


(B7) 


The weight flow out of the bearing periphery is obtained from the 
solution of a simplified Navier-Stokes equation. If the inertia terms 
are assumed to contribute little to the pressure gradient compared to 
the viscous force terms, the equation in cylindrical coordinates (fig. 1) 
is 



(B8) 


If p Is assumed constant in the z-direction, equation (B8) can be 
integrated twice to give 


= - z 2 + Bz + C 
H dr 


where the constants B and C are determined by the conditions that 
u = 0 for z = 0 and for z = h, so that 


u 


.££><*- h ), 
(i dr 


(B9) 


The weight flow w 2 is given by 

w 2 = 3600(24jtrp) /: u ^ (i 11 dr 

or 

dp = _ ^2 

dr 3600(24:trp)h 3 


(BIO) 
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By using equation (B2) and assuming isothermal flow, equation (BIO) can 
he integrated from r s to r, which gives 


P 


2 


-UW 2 RT a r p 

— - In — + pf 

3600itgh 3 r s s 


(Bll) 


With p = p a at r = r a , w 2 is given hy 


w 2 = 


5600 it gh 3 
pRT a in ^ 

1 B 




(B12) 


By using equation (B12), (Bll) can he rewritten as 



or 



(B13) 


Figure 17 is a plot of p/p B against r/r 8 for several values of Pg/Pg 
and r a /r a . A reasonable approximation is obtained hy taking p to vary 

with r as a straight line. With this simplification the expression for 
p becomes 



P - P a = 




(B14) 


The force F-g 


exerted on the upper plate of the hearing is given hy 
yr a 

Fg = 2it / (p - p a )r dr + jcr^(p s - p a ) (B15) 

Jt b 


where the second term on the right is the contribution of the pad pressure 
which is considered to he uniform. Using equation (B14) in (B15) and 
integrating give 




V 


t 


4 


% “ 3 ^ r s + r a r s + r a)(Ps " ^a^ 


(B16) 
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V 


The weight flow W]_ can he related to the weight flow W£ through 
the law of mass conservation: 


w 




3600 


^■(Ps^s) 

at 


or, with V g = nr^(h + d) , 


V 1 - v 2 
3600nr| 


= (h + d) 


d P S 

dt 


dh 
p s dt 


For reversible adiabatic compression. 


Ps Ps,i 



and 




dps 

dt“ 


1 

r 


Ps , i 

p Vr 
*8,1 



<3-Ps 

dt 


(B17) 


(B18) 


Substituting equation (B18) into (B17) and simplifying give 


1 

W 1 ~ w 2 / f \ r fh. + d df dh\ 

3600nr| " Pa ’ i V f J V rf dt d V 


(B19) 


The further assumption is made that the kinetic energy in the orifice flow 
is completely converted to internal energy in the bearing pad so that 
T s,i - Substituting for w-j_ and wg from equations (B6) and (B12) 

and further reducing give 

r 


dh 

dt 


"F^ 


h + d df 

r dt " \s t 


r-i 
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12Aa/r% f r( 
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k 2 r 


12 ^ ( f 2 - J ' 113 


2 , a 

r^u In — 

B r s 


r-i\2 

f r ) ~ 


(B20) 



16 


NACA TN 4095 




Similar reasoning for isothermal compression in the pad leads to 


dh _ 1 

dt " " f 


(h + a.) § - 


12AVe^ i 




(B21) 

Either equation (B20) or (B2l) is used as one of a pair of simultaneous 
differential equations in f and h. Another is obtained by equating the 
product of the mass and acceleration of the bearing to the forces acting 
upon it: 


W d 2 h 
12g ^2 


= (F b - F) 



(B22) 


where a is the ratio of drag to velocity. In the present analysis a 
is always taken to be zero. Substituting from equation (B16) gives the 
second equation in h and f : 

0 ' T + V, + 4)Pr ( f - 5 ) - ¥ * < B23 > 

Equations (B20) and (B23) are the principal equations of the present 
analysis . 
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TABLE I. - BEARING PARAMETERS 
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Clearance parameter 
(h - hfj/Chi - hf) 


HACA IS 4095 


V 



(a) Load variation. 



(b) Bearing clearance. 



Time, t, sec 


(c) Clearance parameter. 

Figure 2. - Load variation and its effect on motion of 
bearing. 



Time, t, sec 


Figure 3. - Clearance variation for basic unstable bearing configuration 
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.0018 
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(a) Bearing clearance. 

Figure 4. - Effect 
at ion of bearing 

of several values of final bearing load on time vari 
clearance and clearance parameter . 
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0 .002 .004 .006 ~ .008 .010 .012 

Time, t, sec 

(b) Clearance parameter. 

Figure 4. - Concluded. Effect of several values of final bearing load 
on time variation of bearing clearance" and clearance parameter. 
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Time, t, sec 

Figure 6. - Effect of variation, of orifice area on clearance 
parameter . 
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Pad radius, r s , in. 

(b) Weight flow. 

Figure 8. - Effect of pad radius on pad pressure and weight 
flow. Fi = 500 pounds j hi = 0.002 inch j r a = 6 inches. 
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Figure 10. - Effect of variation of air temperature on clearance 
parameter . 
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Figure 11. - Effect of variation of mass on clearance parameter. 




Time, t, see 

Figure 12. - Effect on clearance parameter variation due to linear load application over time Interval. Fjl “ 500 pounds; ^ 0,002 
Inch; r a “ 6 Inches; r B ■ 4 inches; A « 0.003 inoh; d * 0.008 Inch; Ff — 550 pounds. 
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Figure 13. - Comparison of effects on clearance parameter of adiabatic and 
isothermal compression in bearing pad. 



hf ) / (hi - h-p ) , dimensionless 
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Figure 14. - Effect of variation of initial loading on clearance parameter. 
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Radius ratio, r/r g 

("b) Bearing radius- ratio r a /r 8 = 1.5. 

Figure 17. - Local pressure ratio as function of radius ratio 
for several conditions . 


NACA - Laniley Field, Va. 







